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We calculate the equation of state of nuclear matter based on the general analysis of the grand 
canonical partition function in the S'-matrix framework. In addition to the low mass stable particles 
and their two-body scattering channels considered earlier, the calculation includes systematically all 
the higher mass particles and their exited states as well as the scattering channels formed by any 
number of these species. We estimate the latter contribution by resonances in all the channels. The 
resulting model-independent virial series for pressure gets substantial contribution from the heavy 
particles and the channels containing them. The series converges for larger values of baryon density 
than found earlier. 
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The equation of state of nuclear matter in the hadronic 
phase is of central interest in the study of supernovae 
and neutron star matter in astrophysical context and of 
heavy ion collisions in the laboratory. There are micro- 
scopic frameworks to deal with the nuclear strong inter- 
action involved here, such as the variational 0,0 and the 
Brueckner-Goldstone methods, using NN poten- 

tials based on experimental scattering data. There are 
also empirical expressions 0| for the Helmholtz free en- 
ergy using a model for nuclei. However, these approaches 
are model-dependent, involving some theoretical uncer- 
tainties. 

A more direct and simple approach based completely 
on observables was suggested long ago by Beth and Uh- 
lenbeck [1, Q, who expressed the second coefficient in 
the virial expansion for pressure of a gas in terms of two- 
body scattering phase-shifts of its constituents. But it 
does not shed much light on the virial series as a whole, 
so that a straightforward calculation up to second order 
may not be trusted, except for very dilute systems. 

A definitive step in this direction was taken by Dashen, 
Ma and Bernstein [l^l, who gave the complete S'-matrix 
formulation of statistical mechanics. Their work is the 
logical completion, though in a highly non-trivial way, 
of the calculation of the second virial coefficient cited 
above. Starting from a few 'elementary' species in the 
definition of the grand partition function, its calculation 
ends up with a multi-species system, with all the bound 
states treated on the same footing as the original 'elemen- 
tary' species. The complete formula involves two sepa- 



rate pieces. The first piece is a sum of ideal gas terms, 
one for each of these stable species. The second piece 
originates from all possible scatterings of any number of 
particles from these species; it is in the form of a virial 
sum over all the scattering channels of terms involving 
their S'-matrix elements. The convergence of the virial 
series is determined by the fugacities and the binding 
energies of these species. 

An attempt to include a bound state heavier than the 
two-nucleon system (deuteron) in the Beth-Uhlenbeck 
scheme has been made by Horowitz and Schwenk [ll|. 
They incorporate the alpha particle explicitly in the def- 
inition of the grand partition function itself along with 
the proton and the neutron and expand it to second or- 
der in their fugacities to get both the ideal gas and the 
scattering terms from these particles to that order. As 
expected, the large binding energy of the alpha particle 
enhances its formation in the system lowering the pres- 
sure. 

In the present work we consider the complete expres- 
sion for the grand partition function [l^, as applied to 
nuclear matter, including, in particular, all the massive 
species and the scattering channels formed by them. The 
only major approximation is that of resonance domina- 
tion of S-matrix elements for scattering channels with 
massive particles. As shown by Dashen and Rajaraman 
it gives rise to additional ideal gas terms from 
these resonances. 

We consider nuclear matter, consisting of protons and 
neutrons as two independent species of 'elementary' par- 
ticles, interacting strongly in the limit of isospin symme- 
try. The object of study is the grand canonical partition 
function for this system. 
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where H is the total Hamiltonian, Np,n the nucleon num- 
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ber operators, (3 the inverse of temperature T of heat 
bath and /ip_„ the nucleon chemical potentials. The 
trace is taken over any complete set of states of all pos- 
sible numbers of nucleons. Denoting the fugacities by 
Qp = e^'^" , = e'^''" , the full trace can be decomposed 
as 

oo 

2= E C^CTrz,Ne-^", (2) 

Z,N=0 

where Trz.N is now taken over states of Z protons and 
N neutrons. For small and C„, hiZ can be expanded 
in a virial double series, 

oo 

lnZ= E MnCpCu- (3) 

Z,N=1 

Our task is to calculate the virial coefficients Az^n- 

We work in natural units, ft = c = 1. For later ref- 
erence, we note here the virial expansion for the ideal 
quantum gas of bosons and of fermions, 

lnZ(°) =^gvj^ ln(lTCe-^'''/^") 

= ^(c±^ + ^ + --] (4) 
A3(m) 25/2 + 35/2 + J ■ 

Here and below the superscript (0) will indicate ideal 
quantum gas. The upper and lower signs correspond re- 
spectively to bosons and fermions. The degeneracy of 
the single particle states is denoted by g V is the 
volume of the system and X{m) = y^2TT/{mT), the ther- 
mal wavelength of a particle of mass m. The integral 
representation above is valid for all values of C, but its 
logarithmic branch point at C = ±1 makes the series con- 
verge only for |C| < 1. 

Though Coulomb interaction is assumed absent in nu- 
clear matter, we shall include its effect on the binding en- 
ergies of nuclei needed below. Our calculation can then 
be readily applied to a physical system, such as neutron 
star matter, by including the effect of added electrons 
necessary to make it electrically neutral. 

We now follow Dashen et al to write the grand 
partition function as the sum of two types of terms. 

In Z = In Z^pll^. + In Z^cat , (5) 

corresponding to contributions from stable, single par- 
ticle states and (multiparticlc) scattering states respec- 
tively. 

Let us first concentrate on the particle piece, InZ^J^^. 
If Z protons and N neutrons form a bound state (nu- 
cleus) of mass number A ^ Z + N, it has mass Am and 
energy. 



in the ground state, where p is its momentum and Bz.n 
the binding energy. From the condition of chemical equi- 
librium among different species, its chemical potential is 
pLz^N = Z^Lp -\- Nfin- Further, all these nuclei, particu- 
larly the heavy ones, have a large number of dense ex- 
cited states above their ground states, which are stable 
in the absence of electromagnetic interaction. (Actually 
their radiation widths arc very small, of the order of eV.) 
Accordingly we split the particle piece in Eq.(5) into two, 

lnZ^::t, = lnZW+lnZW, (7) 

denoting respectively the contribution of the ground and 
the excited states. 

The first term in Eq.(7) is a sum of ideal gas terms, 
one for each of the ground states of all the nuclei, 

InZW =Tl^E5 [t-^ In (1 t e-^('''/2^"-^- "-^-'")) 
z,N J ^'^'^> 

(8) 

where Z and count the number of protons and neu- 
trons in these species. As in Eq.(4) the =p sign correspond 
to nuclei with A even and odd, obeying Bose and Fermi 
statistics respectively. The sum, of course, includes the 
original elementary particles, namely the proton and the 
neutron with their Bz,n = 0. 

The second term in Eq.(7) gives the ideal gas terms 
for each of the elements in the set of excited states of 
the nuclei. Here the individual gas terms are actually 
suppressed by an additional Boltzmann factor, e^^^ , E 
being the excitation energy above the ground state. But 
the density of excited states is rather high in heavy nuclei 
with ^ > 8, say and can, in fact, more than compensate 
for this suppression. We write the contribution of the ex- 
cited states of a single nucleus as an integral over E of an 
ideal gas term multiplied by their level density uj(A, E), 
which we take as [HI, [l^ , 

Lu{A, E) ^ , (9) 

obtained from the Fermi gas model of non-interacting nu- 
cleons within the nucleus. We adopt the empirical value 
for the parameter, a = A/S (MeT^)~^. We thus get 

/ 

InZj^'J =TVJ29 dEu;{A,E)x 

where the prime on ^ denotes exclusion of the light nu- 
clei {A < 8) from the sum. The lower limit Eq is de- 
termined by the beginning of the excited states and the 
applicability of the level density formula at low energy. 
The upper limit E^ is the smallest of separation energies 
of any particle within the nucleus, beyond which lies the 
continuum. We take Eq = 2 MeV and Eg = 8 MeV. 
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We next consider the scattering piece, \nZscat, which 
is formally given by [To| , 



InZ 



scat ^ ^ 



(11) 

where the sum is over all scattering channels, each having 
its chemical potential ^ and formed by taking any num- 
ber of particles from any of the stable species (proton, 
neutron and nuclei in their ground and excited states) 
and the trace is over all plane wave states for each of 
these channels. S is the usual scattering operator and A, 
the boson symmetrization and fermion antisymmetriza- 
tion operator. The subscript c denotes only the con- 
nected parts of the expression in parenthesis. To restate 
Eq.(ll) explicitly for nuclear matter, we have to charac- 
terize the scattering channels. 

Consider the set of channels, in which the constitut- 
ing particles have a total number of Zt protons and Nt 
neutrons with mass number At = Zt + Nt- Let a denote 
all other labels required to fix a channel within this set. 
Clearly the total mass and the chemical potential is in- 
dependent of fT, depending only on Zt and Nt- If P is 
the total momentum of all the particles in a channel, its 
non-relativistic energy is given by 



E 



Zt,Nt,a 



2Atm 



B 



Zt.Nt.a + e , 



(12) 



where e is the kinetic energy in the cm. frame. Here 
Bzt,Nt,u is the sum of individual binding energies of all 
the particles in the channel. With the channels so char- 
acterized, one can integrate over P to get 



In Zscat^ V 



Zt,Nt 



XHAtm) 



j^de e-^^^Trz„N,,, (^^5-i(e)^5(e)^(19) 

the trace being now restricted to the channel {Zt, Nt,(j). 

Eq.(13) brings out the nature of convergence of the 
virial series. Here the integral depends mildly on tem- 
perature. Then the factor e'^^^* ™* " makes channels with 
large binding energies more important, determining the 
leading behaviour of the series. For an estimate of the 
convergence domain, we take Bzt,Nt,a — b{Zt + Nt) for 
such channels with 6 ~ 8.5 MeV, the binding energy 
per nucleon of moderately large nuclei. We thus see 
that Eq.(13) may be regarded as a series in powers of 
g/3(Atp-i-6) pJ3(M^+b) ^ rather than simply of and (n, 
as appears from the original definition (3). If we ignore 
the excited states of the nuclei, the series converges for 
ICpI J ICn| < e~^^ , which, as expected, is smaller than 
that for the ideal gas. Below we shall find a more realis- 
tic criterion of convergence. 

Having discussed the convergence of the virial series, 
we now try to calculate hiZgcat from all the scattering 
channels. It is convenient to divide the channels into light 



ones, consisting of low mass particles {A < 8, say) and 
heavy ones, containing at least one high mass particle 
{A> 8), so that we write 



lnZ„ 



InZ, 



light 



In Z, 



heavy i 



(14) 



as the sum of contributions from the light and the heavy 
channels. The integrals in Eq.(13) over the S-matrix ele- 
ments appear difficult to evaluate at first sight. But the 
presence of the Boltzmann factor simplifies the task by 
limiting the range of integration to rather low energies, 
if the temperature is small enough. 

Generally speaking, two-particle channels are expected 
to dominate over multiparticle channels with same Zt and 
Nt from the binding energy consideration. As in Ref . [ll| 
we neglect all light channels, except the three two-particle 
channels, namely, NN, Na and aa to get 

In Zught = In Znn + In Z^a + In Zaa (15) 

We shall write below the individual contributions in 
terms of phase shifts in the respective channels. 

The restriction of the integrals in Eq.(13) to the low en- 
ergy region is all the more helpful for estimating the scat- 
tering contribution from the heavy channels. The cross- 
sections in these channels are known experimentally to 
be dominated by a multitude of narrow resonances near 
their thresholds. The S'-matrix elements may thus be 
well approximated by these resonances. We then have 
the elegant result by Dashen et al [H, [l^ : The corre- 
sponding partition function becomes that of a number 
of ideal gases, one for each of these resonances. Assum- 
ing their level density to be of the same form as for the 
excited states [l^, we thus see that \nZheavy is again 
given by the same formula (10) as for InZi?, with the 
_B-integral now running from to Er, the end of the 
resonance region. We take Er ~ 12 MeV. We also termi- 
nate the series over Zt and Nt at the same values that we 
choose to do for Z and N in Eqs.(8) and (10). Then the 
two contributions can be combined together by extending 
the variable E from Eq to Er, 

lnZi°) + lnZ,,ea„j; = r.h.s. of Eq.(lO) with E, E,. . 

(16) 

Let us summarize our result at this stage for the grand 
partition function of nuclear matter as, 

In Z = In Z^^,) + (In Z^^J + In Z^eavy) + hi Zi.ght (17) 

with the terms given by Eqs.(8), (16) and (15) respec- 
tively. It is observed that the first two terms represent 
simply the ideal gas terms. 

We now write explicitly the contributions from the 
three light, two-body channels retained in Eq.(15) for 
In Ziight ■ These may be recovered from Ref. [ll[ (see [Ol , 
however), but the master formula (13) gives them imme- 
diately. As we are considering elastic two-body scatter- 
ing, the trace in Eq.(13) becomes a sum over the deriva- 
tive of phase-shifts of the appropriate partial waves. It 
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gives formulae of the same form as derived by Beth and 
Uhlenbeck 0] for the second virial coeiEcient for particles 
without spin and isospin. The results for each of these 
channels can be expressed in terms of an integral of the 
form, 

1 /"^ 

Aab(/3) = — de e-P' ^^^(e), (18) 

where A and B are the scattering particles and 5ab is 
the sum over phase-shifts of the relevant partial waves. 
We shall also indicate the isospin / for the NN channel 
by a superscript on A and 5ab- 
For the AW channel, 

InZjvw = \3{2m) ^^^'^ ^ CI)^''nn 

+CpCn(-3 + A^li + Ay)} (19) 



with 



V(e)= ^(2J+1)<525+1^/^ 



(20) 



S,L,J 



Here the contributing partial waves are determined by 
/ through the requirement of antisymmetry on the total 
wave function of the AW system. Thus 

^^NN = ^Iq + ^in + 3(^3 n + 5(^3 + 5(5i „ H 

•JO ^0 A t^2 J-J2 

= ihg^ + ihp^ + 3(53^^ + 5(53^^ + 7<53 ' 

(21) 



£'3 



The term with —3 in Eq.(19) arises from partial integra- 
tion in Eq.(13), as the phase-shift for the ^Si wave (con- 
taining the deuteron bound statcj_at threshold is not 
but TT, by the Levinson theorem [l8| . 
For the Na channel [l7j . 



InZwa 



V 



where Ba is the binding energy of the alpha-particle and 



'5jva(e) = ^(2J+l)(5L,(e) 



L,.J 



2Ss,,, + 2Sp,., + 4(5P3 + 4^3 + 6Sd, 



/2 



(23) 



Lastly we consider the aa channel with both spin and 
isospin equal to zero, where 



lnZ„ 



V 



A3 (8771) 



with 



<5«a(e) = ^(2L + 1) 5L{e) = Sg + 5Sd + 96c 



(24) 



(25) 



The integrals Aab{P) have also been evaluated in 
Ref. [m with the available phase shifts in the low energy 
region for all the scattering channels. 

Our result (17) for the grand partition function is not 
in the form of the virial series (13). However, we may 
expand each of the ideal gas terms in it again in a virial 
series following Eq.(4) to get the first two terms as 



Z,N 

InZi? -f InZ 



9 



heavy 



Cz,N ± 



25/2 



9 



Z.N 



fl Cz,N ± /2 



Cz,N 
25/2 



(26) 



(27) 



in powers of 'effective fugacitics', Cz,jv 
similar to those already appearing in Eq.(13). The A- 
dependent constants, /„(A), n = 1, 2, • • • 



fn{A) = / dEuj{A,E) e 



-ni3E 



(28) 



are generally large for small 77, if the temperature is not 
too low, but decreases steadily with increasing n, so that 
they make the series (27) converge faster. The ideal gas 
pieces may now be readily calculated, given the binding 
energies of nuclei. We take these binding energies from 
Ref. [101 , where they are available for some 9000 nuclei in 
their ground states up to Z = 135. 

The pressure P and the nucleon densities ob- 
tained from the familiar formulae. 



P = T- 



'V 



d InZ 



i=p,n (29) 



V.T 



at 



and the total baryon density 77b from 77^ = 77p + 
different temperatures and fugacities. 

We now consider symmetric nuclear matter, (p = C" = 
(20| . At each temperature the numerical evaluation of 
density and pressure shows a steep rise in both the quan- 
tities, as we increase the fugacity beyond a certain value. 
We may rely on our evaluation for fugacitics below the 
onset of such rise. It turns out that the allowed range 
of fugacity so determined is also given by half the radius 
of convergence of the virial series obtained earlier in the 
absence of excited states of nuclei, that is, ( < {l/2)e~^^. 
In Fig.l we plot pressure against density for T = 2, 5 and 
10 MeV; where we draw the curves for this range of ( for 
the first two temperatures, while for the third we do so 
for a smaller range to remain well below the saturation 
density. It is found that for fugacities (densities) consid- 
ered here, nuclei heavier than A = 15 do not play any 
significant role. 

For comparison we show in Fig.l the results of 
Horowitz and Schwenk (llj . who calculate Z including 
only one heavy species, namely the alpha particle, besides 
the original proton and neutron and the low mass scat- 
tering channels formed by them, denoted here by Zught- 
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FIG. 1: Plot of pressure of symmetric nuclear matter against 
density at T = 2, 5 and 10 MeV. The solid curves represent 
our calculation while the dashed ones are from the virial cal- 
culation of Ref.flJ], both curves drawn for the same range of 
fugacity at each temperature. The dotted curve with crosses 
(data points) show the microscopic calculation of Ref.[l[ for 
T = 10 MeV. 

It is observed that our results for pressure, besides being 
somewhat lower at higher densities, extend over a wider 
region in density than those in Ref. [ll| . The difFcrcnce is 
to be attributed to nuclei heavier than the alpha particle. 
Also shown are the results of the microscopic calculation 
of Ref.[l| up to a density, beyond which the pressure 
turns out to be unphysical up to a significant region, the 
compressibility there becoming negative. 



Our numerical evaluation shows that the contribution 
from light scattering channels is about an order of magni- 
tude smaller than that from the ideal gas terms. We are 
thus led, to a good approximation, to a modified statis- 
tical equilibrium model for nuclear matter at moderate 
temperature, where we include the contributions from 
not only the nuclei and their excited states as in the 
conventional model, but also the heavy scattering chan- 
nels, which reduce again to ideal gas terms, once the 
^-matrices in these channels are approximated by reso- 
nances present in the low energy region. 

It remains to calculate other interesting quantities, like 
the relative abundance of nuclei and the energy density 
and also consider asymmetric systems. Further, the re- 
sults may be extended to higher temperatures by includ- 
ing pions in the initial set of species. It is, of course, 
known that chiral symmetry restricts the pions them- 
selves to contribute significantly to the pressure [2^ . But 
the contribution from their scattering with the nucleons 
and the a-particles may not be negligible. 

To conclude, we have started from the complete virial 
series for the grand partition function of nuclear matter, 
expressed entirely in terms of observables. Besides the 
binding energies of the nuclei and the density of their 
excited states, the observables include the S'-matrix el- 
ements for scattering of any number of stable nuclear 
particles. While the scattering amplitudes in a few low 
mass, two-body channels are estimated from the experi- 
mental scattering data, those in heavy mass channels are 
assumed to be dominated by low energy resonances. We 
thus derive in a simple way the equation of state of nu- 
clear matter at moderate temperature, that includes all 
significant contributions and, at the same time, is free 
from any serious theoretical uncertainty. Further, the 
results arc valid up to a higher density than found pre- 
viously. 
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